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Abstract 

We rephrase the derivation of black hole radiation so as to take into account, at 
the level of transition amplitudes, the change of the geometry induced by the emission 
process. This enlarged description reveals that the dynamical variables which govern 
the emission are the horizon area and its conjugate time variable. Their conjugation 
is established through the boundary term at the horizon which must be added to the 
canonical action of general relativity in order to obtain a well defined action principle 
when the area varies. These coordinates have already been used by Teitelboim and 
collaborators to compute the partition function of a black hole. We use them to show 
that the probability to emit a particle is given by e~ AA//4 where AA is the decrease in 
horizon area induced by the emission. This expression improves Hawking result which is 
governed by a temperature (given by the surface gravity) in that the specific heat of the 
black hole is no longer neglected. The present derivation of quantum black hole radiation 
is based on the same principles which are used to derive the first law of classical black 
hole thermodynamics. Moreover it also applies to quantum processes associated with 
cosmological or acceleration horizons. These two results indicate that not only black 
holes but all event horizons possess an entropy which governs processes according to 
quantum statistical thermodynamics. 
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1 Introduction 



There are two possible approaches to the gravitational back reaction to Hawking ra- 
diation. The first is to develop a microscopic theory of quantum gravity and to use 
it to calculate the properties of black holes. This program has been partially realized 
in the context of super-string theory Q. The second is to use Hawking's calculation |2|] 
as a starting point to compute the gravitational corrections to black hole evaporation. 
Hopefully these two approaches should meet in some middle ground. 

In the second approach, the back reaction has been addressed along two complemen- 
tary lines. The first is the "semi-classical" theory wherein one first computes the mean 
(quantum average) energy-momentum tensor of the quantized fields propagating on the 
background geometry and then solves the Einstein equations driven by this mean value. 
The solution describes an evaporating geometry characterized by the shrinking of the 
horizon area|| |J In this treatment, the metric remains classical and unaffected by 
the quantum fluctuations of the energy-momentum tensor. 

The second line of attack is to take into account the dynamics of gravity at the level 
of transition amplitudes before performing the quantum average over the configurations 
of the radiation field. Since Hawking radiation is derived from quantum field theory, 
the most natural procedure would be to use Feynman rules. However, the ill-defined 
ultraviolet behavior of quantum gravity has so far prevented this approach from being 
successfully followed. Thus, one has to resort to less ambitious schemes. 

Some interesting insights have been obtained by taking the gravitational back reac- 
tion into account at the classical level, before computing transition amplitudes. To this 
end, the matter action S ma u e r in a given geometry is replaced by the action of matter 
plus gravity S ma tter+gravity A concrete model has been developed by Keski-Vakkuri, 
Kraus and Wilczek (KKW) @, 0, g. It describes the propagation of a massless spher- 
ically symmetric self-gravitating shell in a black hole geometry. Having computed the 
new action, they postulate that the wave function of a shell is given by the WKB form 
e iS m+g _ Tj s j n g these waves, they derive the black hole emission amplitudes of uninteract- 
ing (dilute gas approximation) shells. A similar approach has been used in a Euclidean 
framework in |J following the techniques developed in ]IU| . In this case, the probability 



for the black hole to emit a particle is expressed in terms of the action S m+g of a self 
gravitating instanton. 

The striking result of these works is that the probability for a black hole of mass M 
to emit a particle of energy A is given by 



M^M-X 



N{\, M)e" AA(A ' M)/4 (1) 



where AA(X, M) = A(M) — A(M — A) is the decrease of the area of the black hole 
horizon. N is a phase space (also called grey body) factor which cannot be calculated 
in these approximation schemes. 

Eq. ([[]) replaces Hawking's result Pm-^m-x = Ne~ 2nX ^ K which is characterized by 
a temperature Tjj = k,/2tt defined by the surface gravity n. To first order in A, the 



first law|TT| of classical black hole thermodynamics dE = -^dA guarantees that the two 



expressions for Pm^m-x coincide. This suggests that Hawking radiation and the first 
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law both stem from the same principle. This is far from obvious since the first law 
relates neighboring classical solutions of general relativity whereas black hole radiance 
is derived from QFT in a given geometry. 

The aim of the present work is to reveal their common origin and to establish the 
universal validity of eq. (fl]) . We shall show that both follow from the use of the complete 
action S m+g in place of the matter action in a given geometry. The reason is that the 
emission of a matter quantum can now be viewed as the transition between two neigh- 
boring black hole states. In this transition, the energy is transfered from the hole to the 
radiation, a feature absent in the derivation of black hole radiance based on Bogoliubov 
coefficients evaluated in a given geometry. Moreover, the relevant dynamical quantity 
governing this transition is the difference of actions: S m+g (final) — S m+g (initial) . Then, 
as in deriving the first law[12|, this difference can be reexpressed as a difference of bound- 
ary terms at the horizon. From this, it is easy to show that the dynamical quantity which 
governs emission rates is equal to the change in horizon area divided by 4. 

To demonstrate this result, we make use of the appropriate coordinate system to 
describe processes inducing changes in area. These are the boost parameter G and 
its conjugate variable, the area of the horizon A/8tt. These variables where used by 
Teitelboim and collaborators 1TJ3, O, III], ITS], 113] to compute the partition function of the 



black hole a la Gibbons-Hawking||17|| starting from the appropriate classical action. In 
this case, they showed that one must add to the canonical action (pdq — Hdt) a boundary 
term equal to QA/8tt. In our context, it is through this boundary term that we shall 
obtain eq. (|]) and establish the relationship with the first law. 

We have organized this paper as follows. After a brief review of Hawking's derivation 
formulated in a fixed geometry, we introduce the reader to the variables G and A and 
to their role in boundary terms at the horizon. In particular we emphasize that the 
particular form of the boundary term is dictated by the physical process considered. For 
the emission of particles by a black hole, the appropriate action is the canonical action 
S m+g supplemented by QA/8n. 

This action is used in section § to compute the transition rates of a detector at a fixed 
radius of the black hole. The usefulness of introducing a detector is that the transition 
amplitudes are expressed in terms of the overlap of the initial and final states of the 
detector+black hole complex. Both are stationary eigenvectors of id^ with eigenvalue 
equal to the area of the black hole. In this setting the ADM mass is fixed and therefore 
the time at infinity cannot be used to parameterize the evolution JTB|. Instead one must 
use G time. It is then straightforward to show that the ratio of transition rates of the 
detector is given by e _AA//4 . 

Even though the introduction of a detector is useful, it is not intrinsic to black hole 
radiation. This is why in section |5| we reconsider the KKW model which only makes 
appeal to the action of the particles emitted by the black hole. By using G and A 
we shall recover KKWs result in very simple terms and make contact with the former 
derivation. We shall also show that eq. (0) follows directly from the universal form of 
out-going trajectories in the near horizon geometry and the specification that the field 
configurations be in (Unruh) vacuum. Thus eq. (|I|) applies to all emission processes in 
the presence of horizons, including charged and rotating black hole, cosmological and 
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acceleration horizons (in the latter case the horizon area is infinite, but differences are 
finite and well-defined, see [|10], || |ITfl ). 

The universal validity of this derivation is the main result of this paper. It proves that 
the area of all event horizons determines the gravitational statistical entropy available \ffi 



to these quantum processes. In this we give statistical foundations |^TJ to the relation 
between horizon thermodynamics and Einstein's equations exhibited by Jacobson |22j| . 



2 Hawking Radiation 

In this section we fix the notations and review two standard derivations of black hole 
radiation. In the first derivation, following Unruh[^], we introduce a two level atom 
coupled to the radiation field and whose position is fixed. Then, using Einstein's argu- 
ment, one determines the distribution of massless quanta from its transition rates. In 
this way one only uses basic quantum mechanical rules. The second approach is more in- 
trinsic and closer to Hawking's derivation^]. Black hole radiation is established through 
the Bogoljubov transformation relating in-modes which determine the state of the radi- 
ation field before the collapse (which shall be taken for simplicity to be in- vacuum) and 
out-modes which define the particles emerging from the hole and found at infinity. 
The metric of a Schwarzschild black hole is 

^ = _(i_^ +(1 _^)-V + r W (2) 

We introduce the light like coordinates u and v: 

v,u = t±r* , r* = r + 2M\n\r -2M\ (3) 
and the Kruskal coordinates Uk and Vk- 

U K = --e~ KU , V K = -e m (4) 

K K 

where k = 1/4M is the surface gravity. 

For a black hole formed by the collapse of a spherically symmetric star, the outgoing 
modes, solutions of the Dalembertian equation have the following form near the horizon, 
i.e. for r - 2M < 2M: 

e-^ Y lm (Q) 

<Pujlm = n • (Oj 

\/4tcu r 

Further away from the horizon this expression is no longer exact because of the poten- 
tial barrier which surrounds the black hole. For simplicity, throughout the article, we 
shall neglect the transmission coefficients (also called grey body factors) induced by this 
barrier since they cancel out from the ratio of the transition rates which determines the 
equilibrium distribution around an eternal black hole, see eq. (|T0|) below. This cancel- 
lation also applies to charged and rotating black holes when one takes into account the 
"work term" exerted by the electric potential or the angular momentum of the hole, see 
PI E25|. 
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The modes eq. ([|) are associated with the out-going particles as seen by infalling 
observers. Indeed, for r — 2M <C 2M, the proper time lapse of infalling observers is 
proportional to AUk- This is simply seen by re-expressing the Schwarzschild metric in 
Ur,Vk coordinates: ds 2 ~ —dUKdVx + r 2 d 2 Q. Hence near the horizon, for infalling 
observers, the modes eq. (|5|) have positive frequency. 

Using these modes, the field operator can be decomposed as: 

$ = Yl a ^im^uim + h.c. + ingoing modes . (6) 

By definition, the in- vacuum state, denoted \0u), is annihilated by the a^i m operators. 
In the literature it is often called the Unruh vacuum. The fact that it is annihilated by 
all a^im guarantees that infalling observers experience vacuum conditions as they cross 
the future horizon. 

Consider now a particle detector at fixed radius R of the black hole and angular 
position Q. It has two levels \e) and \g) of energy E e and E g with AE = E e — E g > 0. 
For simplicity we shall take R to be very large (R 3> 2M). Then t is the proper time 
of the detector and AE is the energy gap as measured from r = oo. For smaller r, one 
should take into account the gravitational red shift. 

In the interaction representation, the coupling of the detector to the field $ is given 

by 

H int = rf>(t,R,G) (e +lAEt \e)(g\ + h.c.) (7) 

where 7 is the coupling constant. 

When the detector is initially in its ground state, the state of the system (detector 
plus radiation field $) is \0u)\g)- In the interacting picture, at late times and to first 
order in 7, this state becomes: 

0>g) = \0u)\g)-ilJdte iAEt ^(R,n,t)\0u)\e) 

= \0u)\9) -ilY. I dte^e-^^^alMle) 

= \0u)\g) + ^2B m i,mai lm \0u)\e) (8) 

u)lm 

where C = (l/«) exp^i? + \ ln(i? - 2M)). 

Similarly if the detector was initially in its excited state, the final state would have 
been 

Um J R^/Anuj 
= \0u)\e) + Y, B ~irt,«daJ<>u)\g) ■ (9) 

uilm, 

Since we are interested in determining the population of quanta seen by the detector, 
we only need to compute the ratio of the transition rates. By replacing t by t+m/ ' k in the 



5 



integral governing the transition amplitude B e _» fl)W( j iTO , one obtains B g ^ e = B*_^ g e wAE / K 



for all l,m,u, see |26], |27| . Thus the ratio of the transition probabilities is: 

I R I 2 



15, 
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This corresponds to the rates in a thermal bath at temperature k/2tt = 1/8ttM. Then, 
using Einstein's argument, one obtains that the quanta of the radiation field are dis- 
tributed according to the Planck distribution. 

This derivation is equivalent to calculating the Bogoljubov transformation between 
Unruh modes, eq. (||), and the out-modes defined below in eq. fllTD . Indeed the transi- 
tion amplitudes B g ^ euJ ,i=o and B e ^ guJt i =0 are proportional to the Bogoliubov coefficients 



a u \ and f3^\ when AE = A, see [£3| Pq| - This relationship provides a physical inter- 
pretation of Bogoljubov coefficients as transition amplitudes. A concept that we shall 
rediscuss below. 

We now turn to the second derivation in which the spectrum of emitted particles 
is determined only in terms of solutions of the Dalembertian with a given frequency 
id t = A. There are two positive norm modes for each value of A which identically vanish 
either inside or outside the horizon: 

p— iXu 

V9 A ,+ = ___0(r-2M) , A>0 , 

<fix,- = 4=-*(2M-r) , A>0 . (11) 

Once more we have neglected the potential barrier outside the black hole and for sim- 
plicity we have considered only s- waves (1 = 0). 

The modes (p + define the out-quanta, i.e. those used by a static observer around 
the black hole to describe the presence or absence of out-going particles. We can again 
decompose the field operator into <p± modes: 

$ = a out,\±<P\,± + h.c. + ingoing modes . (12) 

A,± 

By definition the a outy \± annihilate the out vacuum, |0 OM t). 

It is now appropriate to introduce a third set of modes <f>x,± which possess the fol- 
lowing properties. They are eigenmodes of id t = A and are composed only of positive 
frequency modes, eq. (|5]), which define Unruh vacuum. The simplest way to imple- 
ment this last condition is to express the out-modes in Kruskal coordinates: (f\± ~ 
{^■Uk^^^O^Uk)- Since AUk is proportional to the proper time of an infalling ob- 
server, 0a must be the linear combination of tp\± which is analytic and bounded in the 
lower half of the complex Uk plane. Upon requiring also that the 0a± have unit norm, 
one obtains 

^= Vl _Uv. (^ + '" V VU) , A>0 , 
^- = V1-Uv. K+ e " V >; + ) . A>0 . (13) 
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One can verify by evaluating the overlap of and (p^ that the <p\ are linear combinations 
of the 0^ with no 0* component, see e.g. |27| . One can then decompose the field operator 
in terms of these new modes: 

$ = oa±0a,± + h.c. + ingoing modes . (14) 

Then, Unruh vacuum is annihilated by the a\± operators. 

The weights in eq. ( [TB| ) define the Bogoljubov coefficients ot\ and (3\. Their ratio 
satisfies 

\Pm _ -2-kX/k 



\ac\ 

12 



2 



(15) 



Since \(3\\ determines the mean number of out quanta of energy A in Unruh vacuum, eq. 
(li~5|) implies that Unruh vacuum is a thermal distribution of out-particles at temperature 
k,/2tt, in agreement with eq. fllPp. 

We will find it convenient below to use the same argument, but rephrased in coordi- 
nate systems which are regular on the future horizon and which lead to a static metric. 
This second condition implies that the time parameter is proportional to t at fixed r. An 
example is given by the Eddington-Finkelstein coordinates v, r, Q in which the metric 
has the form: 

ds 2 = -(1 -2M/r)dv 2 + 2dvdr + r 2 dtt 2 . (16) 

Near the horizon the metric takes the simple form ds 2 ~ 2dvdr + r 2 dQ 2 which shows that 
v, r are inertial coordinates. Moreover since —dr is proportional to dllx, the momentum 
p r plays the role of the frequency u of eq. (||) and one can translate the analytical con- 
dition implementing Unruh vacuum in terms of r: to obtain a mode, one analytically 
continues ip + in the upper half complex r plane to define its value for r < 2M. As shown 
in [28 1, one immediately obtains eq. (0). In Section 5, it is through the analytical 
behavior in r of the modified modes that we shall determine the corrections to eq. ([15]) . 

In conclusion of this Section, we wish to emphasize the following point. In the 
transitions described in eq. (p]), there is a transfer of energy from the radiation field 
to the detector but the black hole mass stays constant. Similarly, upon computing the 
Bogoliubov coefficients in eq. flTB|), the geometry is unaffected. In these descriptions of 
black hole radiation, there is no transfer of energy from the hole into radiation. The 
notion of black hole evaporation only arises when the mean energy momentum tensor 
of the radiated quanta is put on the r.h.s of Einstein's equations wherein it drives the 
shrinking of the horizon area||-[|5|. 

This artificial two step procedure results from the original sin: to have decided to work 
in a given geometry. As we shall see in Sections 4 and 5, upon working with a recoiling 
geometry, energy conservation will be taken into account at the level of amplitudes, as in 
the Compton effect. Then, the black hole will act as a conventional reservoir of energy: 
when delivering heat to the external world it loses the corresponding energy. 
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3 Boundary terms in the Einstein-Hilbert Action 



In this section, following [|T^, 13), [14|, we introduce coordinates which have an intrinsic 



geometric interpretation near the horizon. These coordinates are the hyperbolic angle B, 
the transverse coordinates along the horizon x\ and the radial proper distance from the 
horizon p. In terms of these coordinates, the metric near the horizon takes the universal 
form 

ds 2 ~ -p 2 d& 2 + dp 2 + 'jijdx^dxi . (17) 

For simplicity we have written the metric for a spherically symmetric horizon, for the 
general case we refer to |14|]. The area of the horizon is A — J dx]_^\ p= Q. In the case 



of a Schwarzschild on-shell solution, O = raf^ and p ~ y 8M(r — 2M) where is the 
proper time at spatial infinity. 

The universal form of the metric in these coordinates implies that if we use them to 
describe a process near a particular horizon, the description of the same process in the 
vicinity of any other horizon will be identical. The developments presented in sections 
[| and [5] for a Schwarzschild black hole thus also apply to charged and rotating holes [15 



and to cosmological and acceleration horizons. Furthermore because these coordinates 
lead to such a simple form for the metric, physical process occurring near the horizon will 
be most simply described in these coordinates. In this section, this will be illustrated by 
considering the boundary terms at the horizon that arise in the Hamiltonian formulation 
of general relativity. In the next sections we shall see that these coordinates are also 
well adapted to describe particle production near event horizons. 

We start the analysis of the boundary terms in the canonical action for matter and 
gravity 

han = f dt' {ir ab g ab + P q-NH- N'H,} (18) 

where g a b is the spatial metric, n ab its conjugate momentum, q and p the coordinates 
and momentum of matter, N and N l the lapse and shift, and H and Hi the energy and 
momentum constraints. In what follows, we focus on metrics g a t, which have two spatial 
boundaries. We suppose that at one boundary the metric is asymptotically flat, that 
is g a b and 7r ab tend to their value in flat space as the proper distance p ~ r tends to 
infinity, see [^9] for the precise behavior one imposes. We also suppose that at the other 
boundary the metric can be put in the form eq. ( fL7|) , see [14] for the precise conditions 



imposed at this boundary. Thus we are considering the class of metrics which, on shell 
and in vacuum, will correspond to one of the asymptotically flat quadrants of an eternal 
black hole. 

Because these metrics have spatial boundaries, it may be necessary to add boundary 
terms to I can in order for it to be stationary on the solutions of the equations of motion. 
We now review these features, starting with the boundary at infinity. Details of the 
calculations will not be presented. They can be found in many papers, see for instance 
P5| P0fl JIEfl [|I4]] . Let us first emphasize that the following considerations are superflu- 
ous for classical physics, i.e. for the determination of the solutions of the equations of 
motion. However, when considering quantum kernels, transition amplitudes or parti- 
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tion functions, they cannot be ignored since the WKB value of the quantum phase is 
determined by the action. 

To identify the boundary terms, following ||13|| , we compute the total variation of the 
action given in eq. fll8D . By total variation we mean that we are considering the linear 
change of I can due to an arbitrary change of all its variables. We note that this variation 
can be performed around any configuration, i.e., on-shell or off-shell. We nevertheless 
impose that the space time is asymptotically flat. In this case, the first order change 
of I can contains three types of terms. First we have the contribution of the bulk, the 
4-geometry interpolating the initial to the final 3-geometry. As usual, this contribution 
vanishes when the reference configuration is on-shell. Secondly, one has the contributions 
due to the change of the initial and final configuration. These determine (as usual) the 
initial and final momenta of gravity. Finally, one has an additional contribution arising 
from spatial infinity. On shell, it is is equal to 

6I r = n °° = t^SMADM (19) 

where too is the (coordinate invariant) proper time at spatial infinity. It is related to 
the lapse function N(r) by too = Jo dt'N(t' , r = oo). 5Madm is the change in mass at 
infinity which is defined by 5g rr for large r. Eq. (fl9|) shows that I can is extremal on 
shell, i.e. its variation reduces on shell to ^ ab 8g a b\{^u a i, only if SMadm = 0, i.e. only if 
one varies among the sub-class of metrics for which Madm is kept fixed. 

If one wishes to consider the Legendre conjugate ensemble in which the asymptotic 
proper time is fixed, one must work with the action S' = I can — t^M adm Indeed, the 
variation of this new action yields the following asymptotic contribution 

SS' r=00 = -bt^MADM ■ (20) 

This term vanishes when one works among the sub-class of metrics for which is kept 
fixed but SMadm is arbitrary. 

We now turn to the boundary term at the horizon. The analysis proceeds in parallel 
with the preceding one. The form of the boundary term is dictated by the fact that 
one requires that, near the horizon, the lapse and shift vanish (or more precisely that 
at p = 0, the momentum n p p and the derivative of the area of the surfaces of constant 
p, d p A, vanish, see jnj) and hence the metric can be put in the form eq. ([TT|) . Upon 
varying eq. (|1|) one finds a boundary term at the horizon 

S jinner boundary = _Q 5A / %lr _ (21) 



Here 8 A is the change in horizon area and G is the hyperbolic angle defined in eq. (|T] 
It is equal to the limit r — > r horizon of N 2 {r)t/2{r — r hor ). For simplicity of writing, we 
have again considered only spherically symmetric 3-geometries. For the general case, see 



[]14 |. Notice that is a coordinate invariant quantity, exactly like too in the asymptotic 



contribution. Moreover, as pointed out in ||12||, A/8n and G are conjugate variables 



exactly like Madm and too. Eq- (23) shows that I can is extremal on-shell for the sub- 
class of metrics which have fixed horizon area. 
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When one wishes to work with a fixed opening angle O, one must work with the 
action S = I can + QA/8tc since its variation yields 

pinner boundary = MQ _ ^ 

S is thus extremal on-shell for the class of metrics which have fixed 0. This action 



must be used when working in the Euclidean continuation of the black hole [16]. Indeed 
regularity of the Euclidean manifold at the horizon imposes a fixed Euclidean angle 
given by 0# = 2tt. The surface term in the action then contributes a term A/4 to the 
partition function which is interpreted as the entropy of the black hole. One of the main 
advantage of introducing this boundary term is to clarify the derivation of this partition 
function which was first considered by Gibbons and Hawking[17[. 



In brief, upon considering dynamical processes occurring around a black hole, there 
are a priori 4 actions which can be considered according to which quantities are fixed 
in the variational principle (or in the path integral) and therefore according to the 
surface terms. Two however are rather unphysical. Indeed fixing both the ADM mass 
and the horizon area is too constraining as can be seen by considering the vacuum 
spherically symmetric solutions for which fixing the ADM mass determines the horizon 
area (in a path integral, the kernel would vanish). For non empty geometries, this 
double specification would impose an unusual non-local constraint on the matter energy 
repartition. Similarly fixing both the time at infinity and the hyperbolic angle is also 
too constraining. Thus one is left with two possibilities: fixing and Madm or fixing 
too and A. A more mathematical reason why these are the only two possibilities (in 
the absence of matter) is that the constraints H = and Hi = viewed as differential 
equations need boundary conditions in order to yield a unique solution and fixing the 



ADM mass or the horizon area but not both provides the required boundary datap4|. 

The choice among the two remaining possibilities is dictated by physical considera- 
tions. If Madm is not fixed but A is, this means that one is considering configurations in 
which the ADM mass is determined by the repartition of matter surrounding the black 
hole while leaving the black hole unchanged. On the other hand, if one fixes Madm 
while letting A to be determined, one is exploring configurations in which energy can be 
redistributed between the black hole and the surrounding matter, but with no change 
of the ADM mass Q. 

Clearly black hole evaporation requires the second situation if one wants to analyze 
what happens at finite r. Therefore in the next section, we shall fix the mass at infinity 
and follow the evolution in terms of and A. The relevant action in this case is 

S = f dt' [n ah g ab +pq-NH -WH^+Q^A/Zk . (23) 



1 To be complete, we should perhaps point out the following difficulty. The horizon is the 2D 
boundary common to all 3-surfaces including the initial and the final ones. Hence, its area cannot vary 
in time. Therefore in order to allow A to vary, one must consider a kind of regularized version in which 
the boundary of the 3-surfaces is arbitrarily close to the horizon, but does not coincide with it. To our 
knowledge, the precise procedure has yet not been completely worked out. We mention here the recent 
works |32| , |33| in which generalized definitions of horizon have been proposed, mainly to allow for 
evolution. However, in what follows, we shall not need an action governing continuous changes in area. 
Only states corresponding to constant areas will be used. 
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In section [5] we shall also consider the region behind the horizon. Thus the formalism 
developed here will not be directly applicable. Nevertheless, it will still be convenient 
to keep the ADM mass fixed and to describe the process in terms of G and A. 



4 Probability for detector transitions in terms of 
horizon area change 

In this section and the following one we calculate the transition rates governing black hole 
radiation when the matter action in a fixed geometry is replaced by the sum of matter and 
Einstein-Hilbert actions. In this new description, the change of the geometry induced 
by the transition process is taken into account through the extremisation of the total 
action. We start with the description of black hole radiation based on the transitions of 
a static detector. 

Our aim is to show that it is the change in area associated with the quantum jump 
of the detector which determines the transition amplitudes. The detector is assumed 
to be at r = R in one of the asymptotically flat quadrants of an eternal black hole. 
We believe however that our analysis also applies to black holes formed by collapse and 
that the mathematical analysis in the two cases should be identical, see however the 
last footnote. As in Section 2, we compute the transition amplitudes in perturbation 
theory to first order in the coupling 7 of the detector to the field, see eq. (Q). As before 
they are given by the overlap of the three free waves. Here these are the radiation wave 
function and the two stationary states of the detector + black hole complex. These 
two stationary states have two interesting properties. Firstly, their eigenvalue is the 
black hole area rather than the detector's energy. Secondly, since the ADM mass is 
fixed, each eigenvalue is entangled to the corresponding detector state. Thus, in the new 
description, the transitions of the detector lead to quantum jumps from one horizon area 
to the other without smooth (classical) interpolation from one stationary geometry to 
the other. 

We first compute the 0-time dependence of the wave functions associated with the 
two states of the detector when their energy is taken into account in the definition of 
the background geometry. This amounts to evaluate twice the on-shell action, eq. (p3|). 
Since the detector is at r = const, both classical geometries are static. Thus the pq and 
^ ah '9ab terms in the action vanish. Moreover, on-shell, the constraints also vanish. Hence, 
at fixed Mad Mi the only term contributing to S is the surface term at the horizon. This 
term is equal to QA g /8iT or QA e /87i where A g (A e ) is the horizon area when the detector 
is in the ground (excited) state. Thus the time dependence of the free (i.e. 7 = 0) wave 
functions are 



This is certainly correct in a WKB approximation. Moreover, since in the absence 
of interactions with the radiation field one deals with stationary area eigenstates, the 



e ieA g /8w ^ 

e iSA e /8TT ^ 




(24) 
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exponential form is exact. Of course, in an exact quantum treatment, the eigen- areas 
may receive quantum corrections, but this will not affect the exponential behavior. 

We now determine the expression for the outgoing modes which replaces eq. ([|). In 
this section, we shall not take into account the gravitational deformation induced by these 
modes. We shall also neglect the effects that the fluctuations of the geometry [[SB], [T?J 
might have. The validity of both hypothesis will be discussed after we have presented 
the mathematical outcome. In the absence of back-reaction and metric fluctuations, the 
modes must be such that they correspond to excitations of the vacuum state defined 
near the horizon. From the definition of O, its relationship to the inertial light like 
coordinates Uk, Vk (such that ds 2 = —oIUkoIVk + r 2 d 2 fl near the horizon) is of the form 
dUic = Ce~ e dO at fixed r. This universal relation exhibits the exponential Doppler 
shift which is the hallmark of horizons (except for extremal black holes). Therefore, we 
make the hypothesis that the new expression is 

W (9, R) = De iuUK = De lCuje ~ e (25) 

in place of e lCu)e K , see eqs. (|| |9]). D is a constant which plays no role. 

We now assume that, to first order in 7 and as in eqs. (|8|, §), the transition amplitudes 
are given by the "time" integral of the product of the three waves. Up to the same overall 
constant (see the discussion after eq. (||)), they then are given by 

Bg^e+ui = JdQ ^BH+g ^*BH+e C > 

Be^g+u, = fdQ** BH+g *BH+e<l>t ■ (26) 



Using eq. (|24f) and eq. ([250 an d by replacing by + in in either amplitude, one 
obtains 

\n i 2 

\ u g^+u\ _ e (A e -A 9 )/4 _ ^27) 

I B e — tg-{-oj | 

This establishes that (A e — A g )/4, the difference of the horizon areas if the detector is 
excited or not, governs the equilibrium distribution of the detector's states. The new dis- 
tribution clearly corresponds to a micro-canonical distribution since we are considering 
exchanges of energy between the black hole and the detector with no ADM mass change 
at spatial infinity. This confirms the interpretation of A/ A as the statistical entropy of 
the black hole since e A//4 enters in eq. ([27]) as the quantum degeneracy of the initial and 
final black hole states. In this we confirm what has been found in |10| , |19| . § . 

Eq. ( |2TD replaces the canonical expression of eq. ( |T0D which is governed by the energy 
change E e — E g and by Hawking temperature k/2ti. To first order in E e — E g , energy 
conservation and the (static version|irT| of the) first law (i.e. dE^etect. = ~dMh i e = 



—KdA/8Tr) guarantee that the new expression gives back Hawking's result eq. (|10|). 



The correction to this first order approximation is governed by the specific heat of the 
black hole. Thus it is completely negligible for large black holes. Therefore the main 
changes from eq. (|To| ) to eq. ( p7|) are conceptual. First, energy conservation is now 
built in through the use of the extremised total action. Secondly, the thermalization of 
the detector no longer reveals that a thermal flux of photons is emitted by the hole but 
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more fundamentally that the detector is in contact with a reservoir whose entropy is 
A/4. Thirdly, the geometry jumps from one stationary situation to another one without 
smooth (classical) interpolation between them. 

Discussion 

In order to reveal the origin of these qualitative changes and to justify the hypothesis 
we made, it is appropriate to rewrite the transition amplitudes in terms of the complete 
system: black hole, detector and radiation field. 

To quantize the whole system requires to consider the Wheeler-DeWitt constraints 

jjBH+detect.+Q^ _ g _ (28) 

When t, the lapse of proper time at spatial flat infinity is fixed, the Wheeler-DeWitt 
equation is supplemented |L2| by the following equation: 

id t \y) = Madm\^) (29) 

Starting from the WDW equation, this shows that one recuperates the notion of a 
Schrodingerian evolution in terms of the (coordinate-invariant) time when one im- 
poses that the 3-geometries are asymptotically flat. Note also that this equation is the 
quantum version of eq. ( p0[) , thereby justifying the classical analysis of that Section. 
Similarly, in the presence of an inner boundary, when working at fixed opening time 



6, the Wheeler-DeWitt equation is supplemented by another remarkable equation|T2[: 



id e \V) = ■ (30) 



57T 



Notice the different signs of A and M in these two equations. This is because A comes 
from an inner boundary and M from an outer one. This flip of sign can already be seen 
at the classical level by comparing eqs. ( |T9"D and (PH). Indeed, eq. flSTf ) corresponds to 
the quantized version of the Hamilton- Jacobi equation 8qS = A which follows from eq. 
(]22f) . In eq. (|30|) , the operator A is defined from the behavior of the 3-geometries as one 
reaches the horizon. When the local WDW constraints are satisfied, A parametrically 
depends on the matter and gravitational configurations from the horizon till r = oo. 
(A simple example of this dependence is provided by the solutions fl24|) of equation (|30|) 
where the eigenvalue A depends on the state of the detector.) 

In our derivation, we first assume that, in the absence of interactions governed by 
7, the energy of the radiation field (whether or not it is in an excited state) does not 
influence the geometry. This amounts to assume that the wave functions of the initial 
and final free states factorize: 

|^n(e)) = \* BH+g (e))®\0u(e)) , 
|*/i„(e)> = |*BH+ e (e)>®|tt $u (e)> . (31) 

This factorization is equivalent to postulate in a path integral formulation that the 
total action splits as a sum: SBH+det+<s> = Ssu+det + S$ wherein the latter is evaluated 
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in the background defined by SBH+det- As shown in pi| , this approximation means 
that the recoil of gravity due to the energy of the radiation field has been taken to 
account to first order only. The corrections to this linear approximation are governed by 
gravitational interactions of the form // TGT where T is the energy momentum tensor 
of the radiation field and G the Green function of the background degrees of freedom ]3"5 



These interactions contain three parts. The first describes the gravitational interactions 
of the photon and the detector. Neglecting these is probably legitimate if the detector is 
far from the black hole. The second part concerns the gravitational self interactions of 
the photon. These will be studied in the context of the shell model of KKW in the next 
section and will be shown to confirm eq. (0). The third part describes the interactions 
among the quanta present in Unruh vacuum. In a dilute gas approximation these are 
neglected. 

In brief, under the assumption of factorizable wave functions, the time dependence of 
the black hole + detector waves are given by eq. ( p4| ) in virtue of eq. (|30|). And eq. (|25| ) 
follows from usual Unruh boundary condition. This last assumption also means that we 
postulate that the near horizon fluctuations have no significant effect on Unruh vacuum. 



One might indeed fear that these fluctuations [36, p7l would destroy its characterization. 



We are tempted to believe however that this is not the case. First the recent body of 



works initiated by Unruh on acoustic black holes ||38|, p9f pleads in favor of this belief. In 
these works indeed, it was shown that there is an adiabatic decoupling between the low 
energy physics governed by the surface gravity and the high energy regime. This explains 
why modifications of the dispersion relation at high frequencies affect neither the low 
energy properties of Hawking radiation nor the characterization of Unruh vacuum. That 



this also applies to the case of near horizon fluctuations is the subject of current work(40 

The second hypothesis concerns the transition amplitudes. To obtain them, one 
should incorporate the interaction Hamiltonian, eq. (0), in the Wheeler-DeWitt con- 
straint, represented here by eqs. (|28| ) and (|30D . To first order in 7, the modified propa- 
gation is expressed in terms of matrix elements of this Hamiltonian sandwiched by the 
free wave functions, as usual for first order in the Born series. This is also true when 
working with the solutions of the Wheeler-DeWitt equation^. Thus to first order in 7 
and up to an overall factor, the transition amplitude is given by 

B g -* e+ul = -i 7 J dO (*BH+e\(V*„\ [HR,S)\e){d\] \* B H+ g )\0u) 

~ -«7 J dQ e - i ( s BH+e-s BH+g ) x e -iS u _ (32) 

In the second line, we have written the phase factors in terms of the action of the black 
hole + detector system and that of the radiation field. This is to emphasize that only 
differences of actions appear in the integrand. Indeed, in the first factor, the action of 
the black hole alone cancels and one is left with the difference due to the change in 

2 In fact the derivation of eq. ( p2| ) closely follows that of transition amplitudes in quantum cosmology 
as performed in . In both cases the wave functions appearing in matrix elements are WKB solutions 
of the Wheeler-DeWitt equation which govern free propagation. Another useful analogy is provided by 
the Unruh effect |23f| . In that case, the given trajectory (which plays the role of the classical geometry 
in black hole physics) is replaced by WKB waves in order to take into account recoil effects fl42||. 
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the detector's state. Similarly, in the second factor, S u is the change in action of the 
radiation field due to the insertion of the field operator at R, 0. Using the stationary 
character of the states for the first factor and Unruh boundary condition for the second 
lead to eq. fl26|) and hence to our central result eq. (|27|) . 

To obtain further insight about eq. Q3"2|), it is interesting to show how it gives back 
the conventional amplitude B g ^ e+UJ obtained in a given background, see eq. @. To 
recover this expression, it suffices to evaluate the difference in actions appearing in eq. 
(|32|) to first order in the energy change: 

S BH +d et. + * (e) _ S BH +d et. + * {e) = ^^{Q) + {E e - E g )d Eiet S BH+det {®) 

= ooU K (Q) + (E e - Eg)t(G) . (33) 

In deriving this, we have first used the splitting of the total action discussed above and 
then Hamilton- Jacobi equations. Indeed, Uk is conjugate to the frequency u and t to 
E detector- It should be stressed that this recovery of the background field phases in the 
limit of small energy differences is a generic feature [34], |19): Whenever one takes 
into account a neglected heavy degree of freedom (here gravity), describes it by WKB 
waves and expands the resulting expressions for light transitions to first order in the 
light change, one recovers the usual background field expressions in which the heavy 
variable is treated classically. 

Moreover this first order expansion "commutes" with the integration over 6. Thus 
it could equally be carried out after having performed the integration. In this second 
form, one is expanding the exponential governing transition rates, eq. (PT[), to first order 
in E e — Eg. In this one recovers the first law of black hole mechanics. This shows that 
the first law is nothing but the Hamilton- Jacobi equation applied to the Euclidean sector 
in the absence of conical singularity: 

d Edet SE u H tL. = d Ede M/^)\M ADM = ~2n/ K . (34) 

In view of the generic character of these features, it is clear that our analysis also 
applies to charged or rotating holes (see however the remark made after eq. ([5])) and to 
cosmological or acceleration horizons. Indeed all that is required is the evaluation of the 
phase factors entering eq. (|32|). In this expression, the second factor e lSiJ presents no 
difficulty: it always encodes vacuum conditions as one crosses the future horizon. The 
first factor is more delicate since it requires to solve Einstein's equations driven by the 
energy of the detector. However in the case of static Rindler like situations, the on shell 
action gets its contribution only from the surface term at the horizon. Moreover, since 
only differences appear in transition amplitudes, eq. fl32|) also applies to the accelerated 
cases. Indeed even though the area of acceleration horizon might be infinite, the change 
in the horizon geometry due to a finite change in the matter energy distribution is finite 
and well defined, see |TtJ for an explicit computation. In particular, it is local in the 
transverse directions x l ± when the change in the matter energy distribution is localized. 
This leads to finite and well defined changes in on-shell actions for the gravity-detector 
system which furthermore gives back the conventional background field result when 
linearized in AE 'detect.- 
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5 The KKW model 

In the previous calculation the change in area was due to the change of the detector state. 
However the existence of a detector is not intrinsic to black hole radiation: the detector 
was only used to reveal the existence of the quanta of the radiation field. Therefore 
we seek for an intrinsic derivation of black hole radiance in which the change in area is 
due to the emission process itself. In this description, the change in area plays the role 
of the energy of the emitted quantum in Hawking's derivation. To this end, we must 
introduce a model for the emitted quanta which takes into account the deformation of 
the gravitational background. The simplest model is that of KKW|| [7], |8| in which one 
makes the hypothesis that in the semi-classical limit and at high frequency the particles 
are described by self gravitating spherically symmetric light-like shells. 

The starting point of the KKW model is the situation analyzed in [53]. It describes 



the entire spherically symmetric space time, solution of Einstein's equations, which re- 
sults from the propagation of a light like shall. By Birkoff's theorem, both outside and 
inside the shell the geometry is Schwarzschild. As in section |] we shall take the outside 
mass Madm to be fixed whereas inside the residual mass M(A) = Madm ~ X depends 
on A, the energy of the emitted shell measured at r = oo. In both geometries the shell 



follows an outgoing light like geodesic, see [43[. In what follows we shall use only the 



inner metric to describe the trajectory, the action and the wave function of the shell. 
This choice will be justified after having presented the results. 

Inside the shell, in Eddington-Finkelstein coordinates the metric is 

dg2 = _ 2M(A) 2 + + r2dQ2 
r 

and the trajectory of the shell satisfies 

dr 



dv 



1 - 2M(X)/r sh 

v(r sh ) = 2r sh + AM{X)\n{r sh -2M{X)) . (36) 

As in the previous section it is a appropriate to introduce a dimensionless time 
parameter defined near the horizon. The light like version of is 

V = K {X)v = 6 + k{\)t + - ln(r - r H {\)) (37) 

where k(X) = 4JV ^ A ^ is the surface gravity and r#(A) = 2M(A) is the final radius of the 
horizon. We shall also use as energy variable the area of the horizon A = 7rr^(A) rather 
then A because A/87T is the conjugate to and to V. In terms of these new variables 
the trajectory is 

V(r ah ;A) = \n(r - r H (A)) + 2k(A)t ~ ln(r - r H (A)) + 0(r - r H (A)) . (38) 

Close to the horizon the log dominates and the trajectory is expressed only in terms 
of quantities locally defined. Therefore this expression characterizes radial trajectories 
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near all event horizons. As an illustration of this universality one easily verifies that for 
a charged non extremal black hole, the trajectory near the horizon also takes the form 
V ~ ln(r — th) wherein the surface gravity does not appear. As in the former Section, 
we make the hypothesis that the near horizon fluctuations do not alter the analytical 
behavior of this trajectory. 

In order to obtain the modes characterizing the quantum propagation of this shell 
we need its action. To obtain the action one could start from scratch, that is from 
the Einstein-Hilbert plus matter action and extremise it. This is the path followed in 
|| |7j. However, since we know the trajectories we can use a short cut to obtain the 
Hamilton- Jacobi action, see also ||. 

Since the inside geometry is static, A is a constant of motion. Therefore the action 
can be written as 

AV 

S(r,V;A) = — + f(r,A) (39) 

where / is the Maupertuis action, that is d r f = p is the momentum of the shell. The 
classical trajectory follows from stationarity of S with respect to A, dj^S = 0. In the 
present case it implies that 8a/ = — V(r; A)/8n where V is given by eq. (j38|). After 
integration from A to A , one obtains 

S(r, V; A) = + [ A ° ^ (ln(r - r H (A)) + k{A)t) - g(r ) . (40) 

A is the area of the horizon in the absence of shell, equal to A = AttM\ dm . With 
this choice, in the absence of the shell the action vanishes. This guarantees that to first 
order in the shell energy A one identically recovers the action of a massless particle in 
the background geometry (at this order it can be taken to be either the inner or outer 
geometry). Notice also that Ao is larger then A since the area of the horizon decreases 
when a particle of positive energy is emitted. 

We have also added an r and V independent integration constant 

r A ° dA / \ 
3( r o) = / ^- (Hro ~ r H (A)) + n(A)r ) . (41) 

At fixed A this function enables the initial momentum p to be fixed arbitrarily and 
determine implicitly the initial radius r or conversely to fix r and determine implicitly 
p . Technically this follows from d ro S = — p (A,r ). 

The action eq. fl4D| ) is well defined on either side of the smeared horizon, that is for 
both r and Tq greater than rn(A ) or both less than ru(A). In each region, we can use it 
to define (in the WKB approximation) the wave function of self gravitating shells. For 
r and r greater than rH(A Q ), one has 

VAA, + = e lS ^ AA ^9(r-r H (A )) . (42) 

This describes the classically allowed propagation of a shell outside the horizon. Similarly 
we can define a wave function 

^ AA _ = e^ v ^e(r H (A)-r) (43) 
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living only in the region r,r < th(A). These definitions are in strict analogy with the 
wave functions (fx± defined in eq. (jTT|) and reduce to them in the limit A —>■ 0. In 
particular the complex conjugation in eq. (|43| ) arises because behind the horizon O time 
runs backwards. Notice that we did not write the prefactors of these waves. To conform 
ourselves to second quantized rules, we should have introduced prefactors such that the 
Wronskian be unity. We shall not pursue this since the prefactors play no role in what 
follows, i.e. the determination of the pair creation probability. 

Between th{Aq) and th{A) one needs a prescription to define the logarithm. As 
explained at the end of section 2, the definition of the action can be used to encode the 
Unruh boundary condition. The analytical specification imposes that the 4>aa,+ mode, 



the analog of 0a,+ of eq. ([131) 



<Ma,+ = e^C^AA) (44) 

be analytical and bounded in the upper half of the complex r plane at fixed V. This 
leads to the globally defined action 

S v (r, V- A, r > r H (A )) = {A ~ A ° )V 

87T 

A> dA 



+ 



A 8tt 



In \r — r H (A) \ + iir9(r H (A) — r) + 0(r — r H (A)) - g(r ) . (45) 



Thus as r goes from r > ru(A ) to r < Th(A) the action acquires an imaginary part 
equal to 

Im Su = ^-^ . (46) 

o 

This simple result is due to the fact that the log in eq. (|4*0| ) comes with an A independent 
weight. As emphasized after eq. ( |38D the origin of this independence follows from the 
universal form of light like outgoing trajectories near a future horizon: V = ln(r— Th{A)). 
Therefore eq. ( |46"D universally follows from this behavior and from Unruh's prescription 
for the analytical behavior of the modes. This is the essential kinematical result of this 
section. 

Using this result, we can write 4>aa,+ as the linear combination 



J A.A, 



VAA, + + e- AA /VAA- ■ (47) 



wherein only the relative weight of <Paa,+ and ^>\a _ has meaning, c.f. the above discus- 
sion about the normalization of the Wronskian. In section |2|, the technique of Bogoljubov 
transformation enabled us to identify the square of the ratio of the weights of <£>+ and 
(p*_ with the probability to emit a pair, see eq. (Jl5|). In the presence of backreaction, 
the same relation still holds for rare and energetic events. Thus we obtain 

Paa = e~ 2ImSu = e" A ^ 4 (48) 

in place of eq. (|T5|) and in agreement with eq. (^) and 0. This result can probably 
be generalized to rotating holes and to cosmological or acceleration horizons. But this 
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will require that one refines the procedure so as to take into account the transversal 
directions. We hope to return to this problem. 

It should be stressed that the justification of the identification of e~ 2ImSu as the 
transition probability is more delicate in the presence of backreaction than in the free 
field theory. This point is discussed below. 

Discussions 

To prepare the discussion, it is appropriate to compare our treatment with the original 
derivation of KKW|| |7|, A first difference with the calculation of KKW is that we use 
A, V instead of A, t as energy and time variables. The passage from one to the other is 
straightforwardly implemented by using the A-dependent Jacobian dA/dX in the action. 
The advantage of using the A, V variables from the start is that it is then manifest that 
the result for the probability of emission P = e~ Aj4//4 is no accident, but follows from 
the universal form of classical trajectories near a horizon. 

A second difference is that KKW work with fixed black hole mass M and with a 
varying ADM mass M + A. In an empty geometry, this can be shown to be mathemati- 
cally equivalent to working as we do with a varying black hole area and fixed ADM mass. 
However as discussed in section ^ the second description reflects better the physics of 
the emission process wherein the black hole loses energy to the radiation while the mass 
at infinity stays constant. Then, when working at fixed ADM mass, the time at infinity 
cannot be used to parameterize the propagation of the shell, see [|18[] for a general proof 
of this super-selection rule. Thus, as in Section 4, it is through the 9 (or V) dependence 
that one recovers the notion of evolution. For non empty geometries, we conjecture that 
it will be mandatory to work with a varying horizon area, i.e. to use the variables A, V 
to parameterize the process. 

Thirdly we have implemented the Unruh boundary condition on the area eigenstates 
0aa,± through the analytical behavior in r as it goes from one side of the horizon to 
the other. In || instead, the Unruh boundary condition is implemented by fixing the 
initial momentum po = —d ro S at V = 0. These two procedures are equivalent. Indeed, 
in both cases, the definition of ln(r — r^) in the action is obtained by imposing positive 
frequency in id r across the horizon. 

The boundary condition of fixed large positive momentum p ^> A also shows the 
relation with the approach of section [|. Indeed the action of a self gravitating shell with 
initial constant momentum po can be shown[] to be 

S(r-V) PQ = [ Po r + h(p )}e- v . (49) 

The assumption of KKW that in the presence of backreaction modes have the form 
_ gi5 P0 j g therefore equivalent to the assumption made in section £|, see eq. (^), that 
the wave function of particles in Unruh vacuum is proportional to e lCue at constant r. 

3 This is obtained by replacing the equation for the trajectory near the horizon V = In ^ 1"^{a) m 

the equation for the momentum p — ^ ^ to obtain p = poe~~ v . Integrating the Hamilton- 

Jacobi equation d r S = p, one finds S Po = pore~ v + F(po, V) where F is an r independent integration 
constant. Then imposing that on the equations of motion the action be stationary with respect to 
variations of po, one finds eq. (pf9| ) with d Po h(po) = rii(ro,Po)- 
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We now turn to the delicate question of identifying twice the imaginary part of the 
action to go from one side of the horizon to the other with the (log of the) probability to 
emit a pair. In addressing this question one faces a double problem. First, due to gravity, 
the emission process is no longer linear. Indeed when two shells of energy A are emitted, 
their action is not twice that of one particle. Therefore, probabilities for multi-particle 
production are no longer obtainable from those for single particle emission as they are 
in a linear field theory. This means that the machinery of Bogoljubov transformations 
no longer applies. The second problem is due to the use of WKB approximations for the 
wavefunctions which are given in terms of the action of a single shell. Because of this 
we are certainly not in a position to describe multi-particle effects, that is higher order 
effects in the tunneling amplitude e~ ImSu . 

We shall now sketch how one can deduce, to leading order in e~ ImSu , the probability 
of particle production from the properties of the wave functions, without resorting to 
Bogoljubov transformations. To present our method, we first return to the analysis in 
the absence of backreaction and consider the following matrix element 

(0u\ J dve iX *$(R = +00, tQjKr, v)\0u) (0e/|<W,a+$M)|0e/) 

= — = <^A+( r > v ) + — <P\-(r, v) (50) 

When r > 2M, it has a simple interpretation. It is the amplitude for an out-particle 
created at (r, v) to be found at X + with energy A. This is a classical process governed 
by the action J r + °° dr'p\(r') where p\ = 2A/(1 — 2M/r) is the classical momentum of the 
particle. The normalization has been chosen so as to describe one out-particle, i.e. to 
have a unit current (Wronskian) for r > 2M. In a linear field theory this normalization 



can be calculated exactly and is given by the factors of a\, see eq. ([13]). It is thus 
proportional to 1 + 0(e _27rA//K ). For self-gravitating shells, this normalization might be 
different but will remain proportional to 1 + 0(e~ 2ImSu ). 

When r < 2M, this matrix element has also a simple interpretation because $(r, v) 
and CLovt^ commute. It defines the amplitude for a pair to be emitted by the black hole. 
Indeed, it defines the transition amplitude from |0jy) to the state $(r, v)a' out x+ \0jj)/a\. 
One member of the pair is outside the horizon and has energy A, the other is behind the 
horizon at (r, v) . /,From this we deduce that the probability P\ to create a pair is given 
by the current f3\/a\ carried by (f)* x+ /a\ for r < 2M. 

This expression for Pa can in turn be used to show how Pa can be expressed in terms 
of the classical action. Recall that the prescription for defining the Unruh mode <p\ + 
for r < 2M is such that its phase, i.e. the action / p\dr, is obtained by analytically 
continuing r — 2M in the upper half complex plane at fixed v. This fixes the ratio of 
the amplitudes 0a+ on each side of the horizon. Thus the probability Pa is given by 

n |0 A , + (r<2M,„)P " 
Note that both the normalization a\ and the relativistic prefactors (47rA)~ 1 ^ 2 cancel in 
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this ratio. This follows from the fact that P\ is given by the ratio of the currents carried 
by 0a,+ on each side of the future horizon. 

Upon taking into account the gravitational backreaction, the above argument still 
applies if two conditions are met. First the WKB approximation of the wave functions 
must be valid. This guarantees that we can use the classical action to evaluate the 
relative amplitude of the Unruh modes on each side of the horizon. Secondly, the wave 
function of the total system, black hole + radiation field, must be (approximatively) 
factorisableQ into the wave function of the shell times the rest. Both require that the 
energy of the shell be large, i.e. X/k^> 1. 

In quantum mechanics such a factorization is a good approximation whenever there is 
a hierarchy in the degrees of freedom: very heavy degrees of freedom that can be treated 
in the WKB approximation (the nuclear degrees of freedom in the case of a Rydberg 
electron orbiting an atom or molecule, the radius of the universe in quantum cosmology, 
the black hole in the present case); moderately energetic degrees of freedom that can also 
be treated in the WKB approximation (the Rydberg electron, heavy matter in cosmology, 
rare energetic particles emitted by the black hole) ; and degrees of freedom that must be 
treated quantum mechanically (the inner electrons, the other matter degrees of freedom 
in cosmology, less energetic particles emitted by the black hole). The very heavy and the 
moderately heavy degrees of freedom propagate semi-classically in the mean potential 
due to the other degrees of freedom (for a detailed treatment of the Born-Oppenheimer 
approximation in a dynamical context, see |0] for atomic physics and |^5[ for cosmology). 
Moreover, the probabilities of the rare quantum transitions of the moderately energetic 
degrees of freedom can be obtained from the sole properties of their wave function. The 
reason is that the wave functions of both the heavier and the lighter degrees cancel 
out from the bra-kets which represent transition amplitudes. An illustration of such 
cancellations is given in eq. (^) where only the difference in actions due to the process 
itself enter in the transition amplitude. 

In brief, when both conditions are met, the former analysis performed in the linear 
theory applies. This guarantees that to first order in e~ ImSu , eq. (|48|) is correct. This is 
also what has been adopted in || [7], |J. 



4 When the shell is very close to the horizon, this factorization probably breaks down because of its 
momentum d r S (= p\(r) in a given background) is arbitrary large (trans-Planckian) since it diverges. 
In the presence of backreaction, the new momentum d r Sjj diverges logarithmicly as r — > ru (A), see eq. 
(fl5|). However, once d r Sjj is much smaller than the Planck mass, the factorization should apply. The 
crucial point is the existence of an intermediate region in which the log is dominating the action Sjj , this 
requires d r Sjj 3> A, and in which the factorization applies, this requires d r Sjj <C Mpi an( ±. In this region, 
the usual analytical characterization of Unruh vacuum still applies [Efl and this is sufficient to obtain 
eq. (48). It is interesting to notice that the same logarithmic behavior (in an intermediate region) 
also explains the absence of modifications to Hawking radiation when one mutilates the dispersion 
relational M 
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6 Conclusion 



Let us further discuss the hypothesis we made. First we supposed that, in the presence of 
backreaction, physics in the neighborhood of a black hole is described by the total action 
Sgravity+matter- Second we supposed that the wave function of the total system factorizes 
into a piece describing the emitted particle and the rest describing the black hole and the 
other matter degrees of freedom. Finally we supposed that the wave functions governing 
heavy degrees of freedom could be approximated by a WKB form. 

The first hypothesis may seem completely evident, but one must recall that there are 
scenarios, such as the brick wall model of 't Hooft, in which it is not necessarily true. 
Furthermore a recent critique of the KKW model does not recover eq. (|I|) because 



its matter model does not derive from an action principle. The validity of the last two 
hypothesis was discussed at the end of the preceding section and in 0. 

Upon making these hypothesis one finds that the probability of the black hole of 
emitting a particle of energy A is given by 



N(\,M)e- AA(x ' M)/4 . (52) 



We have emphasized the universal validity of this result by giving two different deriva- 
tions. Its importance lies in two facts. First it provides statistical foundations to black 
hole thermodynamics since it applies to every quantum emission and since it is gov- 
erned by the induced quantum change (i.e. the recoil) of the horizon area. But, because 
we have not identified the microcanonical states, and because we work in a WKB ap- 
proximation, it is only valid in a mesoscopic sense discussed below. Secondly since the 
derivation of eq. (|52]) also applies to any quantum emission arising in the presence of an 
event horizon, this statistical interpretation is also valid for all event horizons. 



To appreciate the first fact, let us recall that the first law[ll| of black hole thermody- 
namics is a purely classical result which is obtained by comparing two slightly different 
solutions of Eintein's equations. On the other hand Hawking's derivation of black hole 
radiation is obtained by quantizing a matter field in a given geometry. Thus, there is no 
a priori reason why these two results obtained in so different settings should be related 
or even compatible. It turns out however that they are consistent with each other and 
this is the basis for the successful thermodynamics of black holes. 

The origin of this complementarity can be easily understood from the "improved" 
derivation of Hawking radiation presented here. Since energy conservation is imple- 
mented from the start, one takes into account the change of the black hole geometry 
due to the emission itself and one finds that the emission probability is governed by a 
difference of matter +gravity actions S m+g (final) — S m+g (initial) = AA/4. From this 
macrocanonical result one recovers the canonical distribution found by Hawking by ex- 
panding this difference to first order in the energy change A. In the same way the first 
law is also obtained in |12| by this expansion. It is therefore a mathematical necessity 



that they be consistent. 

The other important aspect of eq. ([52] ) is that it indicates that black hole thermo- 
dynamics is not only valid in the mean (in the sense of a large ensemble of processes) 
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but that it applies to every quantum mechanical process. Recall that quantum mechan- 
ics predicts that ratio of the decay rate to the absorption rate of the same quantum is 
proportional to the ratio of the number of final to initial micro-states (at least when the 
interaction Hamiltonian can be treated as a constant). When this is the case, eq. (|52|) 
leads to the identification 

AA/4 = the change of the In of the nbr. of horizon states . (53) 



Moreover, since eq. (|52l) also applies to quantum emissions near all event horizon, eq. 
( |53[ ) also applies to all event horizons. The universality of this identification invites 
three comment sQ. 



Firstly, it leads to the notion of entropy density^, per unit horizon area. Indeed, 
when considering local creation processes rather than global ones described by s-waves, 
the change in area is local [III in the transverse directions x\. This means that a finite 



and localized set of horizon states are affected by the process. 

Secondly, since eq. (^) applies to acceleration horizons, i.e. to near vacuum con- 
figurations having entropy measured by inertial means that can be very small, one is 
lead to conclude that A/ '4 determines the number of states available to the (accelerated) 
system in contact with the event horizon, i.e. with its heat. If this is the case, it implies 
that the black hole area should be viewed as a measure of the entropy accessible to the 
external worldf^, ^j. How to define its 'actual' entropy and whether this latter will still 
be given by A/ 4, as is found for extremal black holes in string theory, [|IJ are important 
questions. 

Thirdly, it should be stressed that the identification expressed by eq. (^) says 
nothing about the microscopic nature of the degrees of freedom which constitute the 
states of event horizons entering in the transition amplitudes we computed. All we can 
say is that when the three hypothesis we made are valid, the number of horizon states 
must obey eq. fl53|). Therefore, eq. (1531) must be conceived as a universal mesoscopic\f)(\, 



21] result, somehow intermediate between the classical laws of thermodynamics and the 
ultimate theory of quantum gravity. 
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